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I.  INTRODUCTION 


By  virtue  of  their  high  sensitivities,  magnetic  field  sensors  employing  SQUIDs 
(Superconductive  Quantum  Interference  Device)  have  been  proposed  for  a  variety  of 
applications  (Refs.  1-4).  They  have  been  used  for  such  diverse  applications  as 
oceanographic  and  geophysical  explorations,  medical  probing  of  the  human  body, 
in  communications,  naval  warfare,  and  for  radioactive  waste  location.  A  particularly  useful 
mode  of  operating  a  SQUID  is  as  a  magnetic  field  gradiometer  (Ref.  5)  which  measures  to 
a  good  approximation  a  spatial  derivative  of  a  component  of  the  magnetic  field.  Thus, 
background  noise  having  long  spatial  correlation  lengths  (such  as  the  earth’s  relatively 
constant  magnetic  field,  for  many  applications)  is  automatically  suppressed. 

If  a  SQUID  gradiometer  can  be  used  to  locate  permanently  magnetic  or  magnetized- 
by-induction  objects  it  should  be  capable  of  seeking,  locating,  and  guiding  a  missile  to  such 
a  target.  It  is  the  type  of  application  that  we  are  interested  in  discussing  here.  Specifically 
we  are  interested  in  evaluating  the  possibility  of  using  the  SQUID  gradiometer  as  part  of  a 
guidance  system  for  an  anti-tank  weapon. 

Figure  1  illustrates  the  concept.  The  superconducting  gradiometers  inside  the 
missiles  sense  the  magnetic  field  gradient  produced  by  the  tank  in  the  earth's  magnetic 
field.  The  seeker  follows  the  gradient  to  the  tank  and  the  warhead  explodes  on  contact. 
The  missile  has  to  be  directed  to  within  about  150  m  from  the  tank,  at  which  point  the 
gradiometer  sensitivity  is  good  enough  to  guide  the  missile.  One  of  the  benefits  of  this 
device  is  that  it  could  be  operated  in  a  shoot-and-forget  mode.  It  could  be  fired  from  a  great 
distance  as  long  as  it  comes  to  within  the  150-m  presumed  distance  or  falls  into  the  150-m 
radius  "bucket"  centered  on  the  tank  when  operated  as  a  mortar.  The  initial  idea  for  this 
application  was  suggested  by  R.  Collins  and  first  described  in  an  IDA  report  published  in 
1988  (Ref.  6). 

The  missile  operated  with  the  gradiometer  is  unaffected  by  common 
countermeasures  affecting  the  operation  of  other  devices  such  as  flares,  camouflage, 
decoys,  etc.  Only  large  chunks  of  metal,  permanent  magnets,  or  magnetic  coils  would 
divert  the  missile.  Thus,  it  is  not  generally  sensitive  to  the  environment. 


Figure  1 


Self-guided  Superconducting  Probe  Weapons  Shown 
Attacking  a  Tank 


m 
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II.  BACKGROUND 


Superconductive  magnetic  gradiometers  have  been  proposed  as  detectors  of 
subsurface  magnetized  materials  both  in  the  ground  and  in  the  water.  One  of  the 
advantages  of  using  a  gradiometer  (as  opposed  to  a  magnetometer)  is  that  the  gradiometer 
automatically  removes  the  large  constant  bias  of  the  geomagnetic  field.  Also,  as  noted 
earlier,  large-scale  spatial  variations  in  the  field  (geomagnetic  noise)  are  suppressed, 
leading  to  potentially  higher  operational  sensitivities  than  are  possible  with  a  magnetometer. 
A  possible  mode  of  operation  for  wide  area  searches  of  sources  of  magnetic  fields, 
magnetic  anomalies,  and  electric  conductivity  anomalies  is  to  place  the  gradiometer  in  an 
airborne  vehicle  and  tow  it  by  an  aircraft.  In  a  passive  mode  such  a  system  could  provide 
large  terrain  coverage  of  the  distribution  of  subterranean  magnetic  anomalies.  In  an  active 
mode,  using  both  a  transntitter  and  sensor,  similar  aircraft-towed  sensor  arrangements  have 
been  employed  extensively  for  subsurface  conductivity  measurements  (electromagnetic 
induction  technique)  (Refs.  7,  8).  A  typical  system  comprises  an  extra  low  frequency 
(ELF)  CW  (200-4000  Hz)  transmitter  on  the  aircraft.  The  EM  fields  penetrating  the  soil  are 
scattered  by  the  spatial  conductivity  inhomogeneities  and  are  detected  by  a  pickup  loop 
housed  in  the  tethered  vehicle. 

If  a  superconducting  gradiometer  could  be  made  out  of  the  high-temperature 
superconducting  material,  it  could  be  employed  as  a  seeker  and  guidance  system  for  the 
destruction  of  tanks  on  the  earth's  surface  or  submarines  under  the  ocean  surface.  Clearly, 
a  guidance  system  gradiometer  sensor  for  both  applications  could  be  made  out  of  a  low 
temperature  superconductive  material.  The  added  complexity  and  weight  due  to  the 
requited  cryogenics,  however,  would  make  the  device  heavier  and  much  more  costly.  For 
the  underwater  application  the  extra  weight  might  be  tolerable,  but  for  the  tank  seeker 
above-ground  application  this  added  complexity  and  weight  would  make  the  concept 
impractical.  The  two  applications  are  similar  conceptually,  although  the  details  are 
different.  In  this  paper  we  will  restrict  the  discussion  to  the  more  demanding  and  difficult 
anti-tank  application. 

Sensitivities  quoted  in  the  literature  achieved  with  SQUID  gradiometers  are  on  the 
order  of  3  x  10~ 1  ^T  m-1Hz-1^  or,  equivalently,  about  10~3  (pT/m)2/Hz  (Ref.  8).  The 
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highest  sensitivity  reported  in  the  literature  of  3  x  10-14  Tm-1  Hz-1''2  was  obtained  by 
Wynn  et  al.  (Ref.  9).  This  intrinsic  noise  level  corresponds  to  the  white  noise  portion  of 
the  spectrum  (f  >  0. 1  Hz);  at  frequencies  below  about  0. 1  Hz  the  noise  power  density  rises 
with  decreasing  frequency,  approximately  in  accordance  with  the  1/f  law.  A  typical 
intrinsic  noise  spectrum  achievable  with  state-of-the-art  SQUID  gradiometer  is  shown  in 
•  Fig.  2. 


Figure  2.  Intrinsic  Noise  Spectrum  of  Magnetic  Field  SQUID  Gradiometer 
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III.  TARGET  DETECTABILITY 


A.  SIGNAL  STRENGTH  AND  SENSITIVITY 

To  determine  the  plausibility  of  this  detection  of  tanks  in  the  earth's  field  we 
estimate  the  size  of  the  expected  perturbation  and  compare  it  to  the  noise  limitations  ^f  the 
SQUID.  Considerations  of  the  earth's  field  are  also  relevant. 

Consider  the  disruption  of  the  ambient  magnetic  field,  assumed  uniform  and  of 
magnitude  B0,  by  a  body  of  magnetic  material.  If  the  material  is  ferromagnetic,  there  is  the 
possibility  that  the  object  will  have  a  permanent  magnetic  moment.  In  general,  therefore, 
the  exact  value  of  the  moment  will  depend  on  the  precise  history  of  the  object  (B-H 
hysteresis)  as  well  as  the  applied  field.  For  simplicity  and  as  a  first  estimate,  for  the  lower 
limit,  the  effect  of  hysteresis  will  be  neglected  and  the  object  characterized  as  a  linear 
paramagnetic  substance  with  permeability  p  »  1 . 

If  the  object  is  further  taken  as  an  ellipsoid  with  the  applied  field  parallel  to  one  of 
the  principal  axes,  then  the  magnetization,  M,  is  also  uniform  and  parallel  to  "he  field.  In 
terms  of  the  demagnetizing  factor  y  for  that  axis,  M  is  given  by  (Ref.  10): 


4kM  = 


(H-DB0 

(i-y)  +  yp 


0) 


For  large  |i,  and  y*  0  ,  the  magnetization  is  essentially  independent  of  the  precise 
value  of  p.  A  long  cylinder  with  the  field  along  its  axis  has  a  value  of  y  =  0.  For  a  sphere, 
y  =  1/3  and  all  axes  are  equivalent.  For  a  sphere  of  radius  R,  *he  total  magnetic  moment, 
m ,  is  simply  the  product  of  the  volume,  47tR3/3,  and  the  magnetization,  M,  or  m  =  BcR3. 

The  total  field  is  the  sum  of  the  applied  local  field  and  the  dipole  field  due  to  the 
induced  magnetization  of  the  "para"  or  "ferro"  magnetic  material.  A  magnetometer  can  be 
used  to  determine  the  field  at  any  point  in  space  and  a  gradient  magnetometer  senses 
changes  in  the  total  field.  The  limits  to  determining  the  field  or  its  spatial  derivatives 
(gradient)  is  set  by  the  noise  in  the  detector  (sensor)  systems.  Instruments  for  measuring 
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the  field  gradient  have  been  developed  over  the  past  decade  and  superconducting  devices 
have  been  reported  to  have  limits  of  the  order  of  10~12  G/cm/Hz1/2.  Thi.  noise  limit  "S" 
may  have  been  due  to  the  technology  used  at  the  time  and  is  used  below  for  estimation 
purposes. 

For  the  spherical  object  described  above,  the  magnitude  of  the  dipole  field  is 


B  =  mjr$  . 


(2) 


where  r  is  the  radial  distance  to  the  dipole.  In  this  crse,  the  gradient 


dB 

dr 


(3) 


This  magnitude  of  the  gradient  must  be  compared  to  the  noise  "S,"  that  is. 


>  S  (Af) 
dr 


1/2 


(4) 


In  terms  of  the  noise  limit,  S,  the  range  of  detection  of  a  spherical  mass  is  given  by 


r  < 


3B0R? 
S  (Af)1/2 


(5) 


An  estimate  of  "r"  for  a  typical  system  is  200  m.  Parameters  used  in  this  estimate 
are  the  mean  field  of  the  earth  and  a  sphere  nf  1-m  radius.  The  mean  field  of  the  earth  is 
approximately  0.5  G,  a  sphere  of  1-m  radius  is  approximately  35  tons  of  iron.  A  system 
bandwidth  of  100  Hz  was  taken. 


Another,  perhaps  ultimate,  limit  is  given  in  terms  of  flux  fluctuations  by 
Sfiux2  =  h  L,  where  h  is  Planck's  constant  and  L  is  the  inductance  of  the  SQUID  loop.  For 


typical  loops  L  =  10-9  h,  a  gradient  limit  produced  by  placing  two  such  loops  a  distance  d 
cm  apart  is  approximately  Sd  =  8  x  lO-14  G/Hz-1/2.  Ford  =  10  cm,  and  a  100-Hz 
bandwidth,  this  provides  an  additional  factor  of  three,  or  rdetecljon  =  600  m.  Clearly,  it  is 
of  great  importance  to  further  define  the  noise  limit  of  such  a  device. 

Use  of  the  dipole  gradient  makes  the  range  relatively  insensitive  to  changes  in  the 
detector.  For  example,  increasing  the  size  of  the  SQUID  loop  increases  the  signal  strength 
proportional  to  the  area  of  the  loop,  but  increases  the  noise  by  the  square  root  of  the  area. 
Thus,  changing  the  area  of  the  loops  by  a  factor  of  100  only  increases  S  by  a  factor  of  10, 
and  the  range  by  approximately  70  percent. 

A  more  severe  limit  than  the  SQUID  noise  may  be  the  fluctuations  in  the  earth's 
field.  This  may  provide  a  more  stringent  limit  than  the  performance  of  the  SQUIDs. 

The  simple  model  offered  above  suggests  several  applications,  two  of  which  may 
be  of  interest  to  DoD  and  are  described  below. 

B  .  MAGNETIC  FIELD  AND  GRADIENT  FROM  A  DIPOLE  SOURCE, 

A  COLLECTION  OF  DIPOLE  SOURCES  AND  AN  EXTENDED 
OBJECT  OR  MULTIPLE  SOURCE 

— *  — >  _ _ ^ 

The  magnetic  field  B(r)  at  point  r  due  to  a  dipole  located  at  the  origin,  as  shown  in 
Fig.  3,  is  given  by: 


3n  (n  •  m)  -  m 


(6) 


where  n  is  the  unit  vector  pointing  in  the  r  direction,  m  is  the  magnetic  dipole  moment. 
For  the  magnetic  dipole  moment  at  the  origin  point  in  the  z  direction  we  have 


BfX)  = 


!m  i  (3n  cos6  -  k) 

rT 


(7) 


where 


where  B0  is  the  static  field  and  K  =  lm  I  is  the  magnetic  dipole  moment  strength.  For 
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spherical  objects  of  radius  Rq  of  uniform  magnetization  K  =  Rq  B0- 

The  components  of  the  field 

B(X)  =  Bxi  +  Byj  +  Bgk  (11) 
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are 


'  lx  B*  =  3Kxz  r'5 
-  =  B,  =  3Kyzr'5. 


I 

The  components  of  the  gradient  dyadic 


> 


) 


where 


Qxy 

Q,,' 

=> 

Q  = 

Qyx 

Qyy 

Qyz 

{Qzx 

Q^y 

QzzJ 

^  Ba  - 


and  the  components  are  given  by 


Qxx  =  3KCZT2  -  5x2z)  r'7 
Qxy  =  -  15Kxyzr'7 
Qxz  =  3K(xr2  -  5xz)  r'7 
Qyx  =  -  1 5Kxyz  r'7 
Qyy  =  SKOr2  -  5y2z)  r  7 
Qyz  =  BKCyr2  -  5yz2)  r'7 
Qzx  =  3K(-5xz2  +  xr2)  r'7 
Qzy  =  3K(-5yz2  +  yr2)  r'7 
Qm  =  3K(2zr2  -  5z3)  r7 


(12) 


(13) 


(14) 


(15) 


For  a  collection  of  N  dipoles  located  at  different  points  r;  and  with  different  dipole 
strengths  k,  we  have 
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<D  =  zB0 


(z  -  z.) 


r  -  r; 


B  =  k 


(16) 


N 

+  3  £  K.(z  -  z.) 
i=l 


(r 


For  an  extended  object,  the  obvious  extension  of  the  above  expressions  is  required. 

C .  LOCATION  OF  MAGNETIC  OBJECT  FROM  G RADIOMETER  SIGNAL  # 

The  purpose  of  this  section  is  to  develop  an  algorithm,  that  can  be  easily 
implemented  on  the  computer,  for  locating  a  magnetic  object  from  a  gradient  of  the 
magnetic  field  at  a  point  remote  from  the  object.  In  subsection  (A)  we  derive  an  algorithm 
locating  a  magnetic  dipole  (object)  in  a  constant  earth’s  magnetic  field.  In  subsection  (B)  • 

we  develop  a  procedure  for  dealing  with  an  extended  object  whose  signature  is  then 
described  in  terms  of  a  multipole  expansion. 

I .  Point  Source  Magnetic  Object  in  a  Constant  Field  0 

Consider  the  magnetic  field  composed  of  a  uniform  field  and  a  dipole  oriented  with 
the  field;  the  magnetostatic  potential  for  the  far  field  is  given  by 


<f)  =  B0  z  -  K  -J  .  (17) 

If  the  field  direction  is  unknown  (or  if  the  dipole  is  not  oriented  along  the  field),  then  the 
only  information  available  is  given  by  the  tensor  gradient  of  the  magnetic  field: 


(18) 
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where  the  M  matrix  is  given  by 


M..  =  cos0  [8..  -  5  r.r.]  +  [rz  +  z  r  ].. 

ij  ij  i  j 


(19) 


or  in  other  coordinates  or  in  the  coordinates  shown  in  Fig.  3. 

M..  =  cosQ  [5.-3  r.r  ]  -  sin0  [r  0  +  0  r  ].. 
ij  1  ij  i  r  1  Jy 


(20) 


where  0  is  the  polar  angle  from  the  z  axis  (the  dipole  axis).  The  prefactor  of  M  in  Eq.  (18) 
involves  the  unknown  distance  and  dipole  strength  of  the  dipole  under  consideration  and 
therefore  (at  least  at  first)  only  the  relative  sizes  of  the  eigenvalues  of  the  matrix  are 
relevant.  [For  a  sphere  of  permeability  |i  and  radius  Rq,  K  =  (p.  -  1)  B0R03/(H  +  3).] 
Note  that  M  has  a  simple  form  if  the  orientation  of  the  dipole  is  known;  this  information  is 
not  always  available  but  the  procedure  below  can  be  simplified  if  it  is  known. 

M  has,  in  general,  three  nondegenerate  eigenvalues  and  eigenvectors.  The 
azimuthal  direction,  <>,  is  always  an  eigenvector,  since  it  is  perpendicular  to  both  the  radial 
and  z  directions  (alternatively  to  the  r  and  0  directions).  Examining  the  M  matrix,  it  is  clear 
that  this  eigenvector  has  eigenvalue 


X.  =  cos0  . 

<p 


(21) 


The  two  remaining  eigenvalues,  corresponding  to  superpositions  of  the  t  and  0  directions 
are  solutions  of 


X  +  cos0  X  -  (1  +  cos20)  =  0. 


(22) 


The  two  values  straddle  the  azimuthal  eigenvalue 


(23) 


i  _  -  cosQ  ±J  4  +  5  cos20 
X±~  2 

[Note  that  since  the  prefactor  is  unknown  the  eigenvalues  of  the  field  gradient  tensor  will  be 
these  values  multiplied  by  3K/I4.]  Except  for  cos  0  =  ±1,  these  eigenvalues  are 
nondegenerate.  By  computing  ratios  of  the  three  eigenvalues  of  the  gradient  tensor,  two 
independent  measures  can  be  taken  of  6.  For  example,  one  can  consider  the  middle 
eigenvalue  divided  by  the  difference  of  the  largest  and  smallest. 

COS0  _  _  middle _ 

J -  1 

V  4  +  5  COS20  oiggest  smallest 


Note  that  the  sum  of  the  three  eigenvalues  must  be  zero  because  the  magnetic  field  is 
divergenceless  (this  is  true  for  an  arbitrary  source  field,  not  just  the  single  dipole  field);  if 
the  measured  values  do  not  satisfy  this  bound,  it  is  a  measure  of  their  error.  The 
eigenvectors  that  correspond  to  the  eigenvalues  of  Eq.  (23)  are  combinations  of  the  radial 
and  0  directions.  Writing  the  (unnormalized)  eigenvectors  as 


v 


± 


(24) 


one  has  the  equations 


“K±  =  -  2  cos0  -  sin0  a± 

X±a±  =  -  sin0  +  cos0  a± 


(25) 


which  determines  a±  in  terms  of  the  eigenvalues  of  the  M,  (or,  equivalently,  having 
determined  the  polar  angle  0  from  the  ratio  of  the  eigenvalues  the  formulas  above,  in  terms 
of  that  angle).  Note  that  having  determined  the  polar  angle,  the  value  of  the  prefactor  3K/I4 
can  be  determined,  which  gives  an  estimate  of  the  range  if  the  rough  size  of  the  dipole  is 
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known.  Having  determined  the  eigenvectors  of  M,  the  radial  direction  can  be  recovered  by 
solving  the  above: 


A 

r  = 


(26) 


The  procedure  for  steering  is  therefore  as  follows: 

*  Compute  the  eigenvalues  of  the  gradient  tensor,  thereby  determining  the  polar 
angle. 

•  Compute  the  eigenvectors  and  form  the  combination  given  in  Eq.  (26). 

•  Estimate  the  range  for  a  typical  size  dipole  via  the  prefactor. 

*  Steer  (obviously  information  from  previous  time  step  is  useful:  Make  a 
Kalman  filter!). 

Unfortunately,  the  sign  of  the  eigenvectors  is  arbitrary  so  that  in  the  absence  of  any  other 
information,  Eq.  (26)  gives  four  possibilities  for  the  radial  direction.  One  can  choose  the 
one  closest  to  the  direction  of  motion,  or  use  putative  information  about  the  location  of  the 
dipole  z-axis  or  location  gravitational  vertical  to  disambiguate. 

2 ,  Extended  Magnetic  Object  (Point  Source  Multipole  Expansion) 

The  analytical  calculation  given  for  the  steering  algorithm  presumes  a  single  point 
dipole  source  for  the  magnetic  field  (in  addition  to  any  constant  background  field).  The 
numerical  studies  of  the  steering  algorithm  assumed  several  additional  point  dipole  sources 
which  could  act  as  decoy  targets.  Whenever  one  of  these  targets  is  clearly  dominant  the 
guidance  algorithm  is  essentially  guaranteed  to  converge  on  the  dominant  point  dipole. 
However,  real  extended  sources  will  have  more  complex  signatures  than  a  simple  dipole. 
As  the  source  is  approached,  higher  multipole  contributions  to  the  field  will  perturb  the 
gradient  tensor  and  at  sufficiently  short  ranges  may  dominate  the  dipole  contribution.  This 
might  lead  at  short  range  to  a  deflected  trajectory.  This  can  be  explored  numerically  with 
extended  sources  or  sets  of  closely  spaced  dipole  sources.  A  complete  analytic  discussion 
is  beyond  the  scope  of  the  present  effort.  An  initial  outline  of  a  multipole  perturbation  of 
the  steering  algorithm  will  be  given  in  this  section. 
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For  an  arbitrary  scalar  potential  O,  the  gradient  tensor  is  given  by 


A  A 

O.  =60 
ij 


A  ^ 

AA  _  A  A 

r  r2 

l—  _J 

+  rr  <t>n  +  <jxt> 

cb. 


r  .  °ee  „  °ecose 


+  r0  + r0j 


O 


q> 

r0  '  2 

r 


aa  aa1 


+  |  4>r  +  r<bj 


sin20  r2  i^sin© 


6  O 

yr6  o 


rsin0  sin0r2 


(27) 


"  0> 

COS0 

'  A  A  A  Al 

4>0  +  0<f)J 

60 

6 

t^sin© 

t^sin2©  . 

The  extended  source  can  be  represented  in  the  "far”  field  by  an  expansion  in  spherical 
harmonics. 


1  P.(cos0)  v  YJO,*) 
0  -  - 1 — -  +  >  a, 


St2 


lm 


ta  3r,+1 


(28) 


where  the  1/3  is  introduced  to  match  the  prefactors  used  earlier. 

One  may  imagine  considering  the  effect  of  a  single  Y im  term  on  the  steering  in  the 
following  way.  The  gradient  tensor  can  be  computed  with  the  dipole  and  perturbational 
term  included.  The  eigenvectors  and  eigenvalues  can  be  computed  and  the  operations  of 
the  steering  algorithm  given  previously  applied  to  the  new  situation.  The  error  in  the 
assignment  can  then  be  computed  as  a  function  of  the  aim  coefficient.  To  first  order  in 
perturbation  theory  these  effects  may  be  summed.  The  algebra  needed  for  this  program  is 
tedious  but  not  difficult  and  will  be  pursued  in  later  work.  For  current  purposes,  only  a 
few  general  remarks  need  be  made. 

First,  the  vanishing  of  the  divergence  of  the  magentic  field  is  equivalent  to  the 
vanishing  of  the  trace  of  the  gradient  tensor  (equivalently,  the  eigenvalues  sum  to  zero). 
This  condition  holds  for  any  multipole  expansion  as  can  be  verified  by  direct  examination 
of  the  equivalent  partial  differential  equation  for  O. 
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Second,  in  the  dipole-only  case  considered  earlier,  the  largest  and  smallest 
eigenvectors  are  contained  in  the  r-0  plane.  Thus,  the  operations  given  above  are 
equivalent  to  determining  the  angle  through  which  the  largest  and  smallest  eigenvectors 
(which  are  orthogonal)  must  be  rotated  to  align  them  with  the  r-0  directions.  For  the 
multipole  expansion  case  with  m  =  0  (no  explicit  <}>  dependence)  this  property  is  preserved. 
The  effect  of  the  higher  order  multipoles  is  to  change  the  required  angle  of  rotation  but  the 
radial  direction  still  lies  in  the  plane  spanned  by  the  largest  and  smallest  eigenvectors.  (For 
infinitesimal  perturbations,  the  relative  size  of  the  eigenvectors  is  unchanged  whenever 
there  are  three  distinct  eigenvalues.)  For  m  *  0,  this  does  not  hold  and  the  true  radial 
direction  will  lie  out  of  the  plane  of  the  largest  and  smallest  eigenvectors. 

The  algebra  for  even  the  simplest  example  of  m  =  0,  1  =  2  (quadrapole)  case  does 
not  appear  to  be  immediately  enlightening  and  will  not  be  reproduced  here.  Progress  in  the 
multipole  expansion  requires  either  a  long  series  of  calculations  to  be  done  by  hand  or  the 
use  of  a  symbolic  manipulation  program  such  as  MACSYMA. 
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IV.  STEERING  CODE  DEVELOPMENT  AND  TESTING 


A.  BACKGROUND 

Using  the  theoretical  basis  laid  down  in  Section  II  we  developed  a  code  which  is 
designed  to  utilize  gradiometer  measurements  to  steer  a  missile.  In  this  section  we  describe 
the  code  and  examine  its  capabilities.  The  code,  which  we  will  refer  to  as  XGRAP  inverts 
the  gradiometer  signal  to  obtain  a  direction  which  points  to  the  source  of  the  field.  It  turns 
out  that  such  an  inversion  is  not  unique  but  produces,  in  general,  four  vectors,  two 
pointing  in  the  directions  of  the  source  or  more  accurately  into  the  hemisphere  in  which  the 
source  is  located  and  the  two  others  point  into  the  opposite  hemisphere.  These  latter  two 
can  be  dispensed  with  easily  on  grounds  that  the  missile  knows  the  difference  between 
forward  and  backward  directions. 

Clearly,  to  develop  a  code  which  can  be  used  in  the  field  to  guide  missiles  requires 
many  steps  including: 

1.  Test  of  mathematical  validity  of  the  code  under  simulated  operational 
conditions,  assuming  pure  signal  (no  noise  or  decoys). 

2.  Determination  of  the  effects  of  noise  both  natural  and  manmade. 

3 .  Determination  of  the  effect  on  results  of  false  targets,  including  decoys. 

4.  Determination  of  the  effect  on  steering  of  ambiguous  signals  (i.e.,  missile 
equidistant  from  two  targets. 

5 .  Determination  of  the  range  of  operation  of  system. 

6.  Simulation  of  all  relevant  physical  processes,  including  flight  characteristics, 
acceleration,  steering  response,  etc. 

In  this  paper  we  touch  on  all  of  the  above  points  with  the  exception  of  (6).  From 
the  examples  to  be  presented  here  it  will  become  clear  that  more  in-depth  analysis  is 
required  to  determine  the  realistic  bounds  of  operability  of  the  proposed  steering  system. 
The  intent  here  is  to  show  feasibility  rather  than  to  attempt  to  conduct  an  exhaustive  study 
of  special  cases. 


B.  PROGRAM  XGRAP-ALGORITHM  STRUCTURE 


The  following  steps  are  taken  in  the  code  to  steer  the  missile  in  a  course  toward  the 
magnetic  source. 

=> 

Step  1.  Calculations  of  the  magnetic  field  gradient  Q  at  the  missile  position  x0,  y0, 
Zq.  This  is  stored  as  a  matrix 


Q(x0,  Y0,  zQ)  — 


fQxx  Qxy  Qxz^l 
Qyx  Qyy  Qyz 
VQzx  Qzy  Qyz  j 


(29) 


Step  2.  Calculations  of  the  eigenvalues  and  eigenfunction  of  Q  xQ,  y0,  Zo  as 
shown  below. 


Eigenvalues 

Order 

Eigenvectors 

X 

+ 

high 

— » 
v+ 

X 

least 

X 

middle 

— > 
v„ 

0 

o 

Step  3.  Calculation  of  cos0  from  equation 

cos6  _  (30) 

4  +  5  cos20 

and  a  check  on  the  eigenvalues  to  see  if 

X+  +  X.  +  Xo  =  0.  (31) 
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Step  4.  Description  of  eigenvectors  in  terms  of  angular  unit  vectors  so  that 


A  £ 

v+  =  r  +  a+  0 

— >  a  £ 

v.  =  r  +  a.  0 


and  a+,  a-  are  determined  from 


X+  =  -  2  cos0  -  sin0  a+ 
X.  =  -  2  cos0  -  sin0  a. 

Two  checks  are  made  at  this  point  to  determine  if 

X±a±  =  -  sin0  +  cos0  a+ 


and 


(32) 


(33) 


(34) 


(35) 


Step  5.  Determination  of  the  direction  of  the  target  Rp  by 


— »  — > 

a  -  v+  -  a+v 

a--a. 


(36) 


Step  6.  The  missile  is  moved  from  point  x0,  y0,  z0  to  point  yi,  yi,  z\  in  the 
— ) 

direction  Rp,  at  which  time  the  process  of  following  steps  1  through  5  is  repeated  and  a 
new  point  X2,  y2>  z  is  reached.  The  process  is  depicted  in  Fig.  4. 


19 


POINT  REACHED  AFTER  SECOND 
STEERING  CORRECTION 


Figure  4.  Illustration  of  the  Steering  Process 


C.  TESTING  OF  CODE  XGRAP 

In  preparation  for  the  steering  algorithm  tests  we  made  a  study  of  the  field  B  and 
=s> 

field  gradient  Q  resulting  from  different  magnetic  dipole  distributions  as  shown  in  Fig.  5. 
The  following  four  figures  (Figs.  A-l  through  Fig.  A-4)  show  some  sample  field  and 
gradient  distributions.  In  these  figures  a  grid  of  10  x  10  x  10  points  is  used  and  the  field 
vector  and  gradient  matrix  were  calculated  at  each  point.  In  each  case  considered,  a  cut  is 
taken  through  the  y-axis  at  y  =  1  and  the  values  of  the  fields  at  the  100  points  in  the  (x,  y  = 
1 .  z)  plane  are  shown. 

The  basic  format  of  the  figures  is  as  follows:  the  top  three  views  in  each  figure  (a, 
b,  c)  show  the  equal  magnetic  field  contour  lines  of  the  components  Bx,  By.  Bz. 
respectively.  The  next  six  views  (d,  e,  f,  g,  h,  i)  show  the  equal  magnetic  field  gradient 
contour  lines  Qxx,  Qyy,  Qzz-  Qxy*  Qxz.  Qyz,  respectively.  The  other  components  are 
redundant  because  of  symmetry.  In  each  figure  the  logarithm  of  the  value  is  actually 
plotted. 


Table  1  shows  the  location  of  the  components  in  Figs.  A-l  through  A-4. 


Table  1.  Location  of  Magnetic  Field  and  Gradient  Components  In 

Figs.  A-1  through  A-4. 


B 

B 

B 

X 

y 

z 

Q 

Q 

Q 

XX 

yy 

zz 

O 

Q 

Q 

xy 

xy 

yz 

Figure  5  shows  the  location  and  relative  strengths  of  the  dipoles  used  in  the 
subsequent  calculations.  The  plane  (x,  y  =  1,  z)  is  also  indicated  by  shading.  In  Fig.  A-l 
the  contour  plots  for  fields  resulting  from  a  single  dipole  of  strength  1 .0  located  at  (0,  0,  0) 
are  shown.  In  Fig.  A-2,  in  addidon  to  the  dipole  at  (0,  0,  0)  of  strength  1.0,  another  dipole 
of  strength  0.1  is  assumed  located  at  (0,  10,  0).  In  Fig.  A-3,  the  second  dipole  is  moved  to 
(10,  0,  0)  and  its  strength  is  increased  to  0.5.  To  calculate  the  results  shown  in  the  last 
figure  of  this  series.  Fig.  A-4,  three  dipoles  are  assumed.  The  first  dipole  of  strength  1.0 
is  at  (0,  0,  0)  the  second  is  at  (0,  10,  0)  and  has  a  strength  of  0.2  and  the  third  is  located  at 
(5,  5,  0)  and  has  a  relauve  dipole  strength  of  0.1. 

It  is  difficult  to  determine  the  effect  on  the  steering  algorithm  from  such  plots.  The 
decision  was  made  to  just  use  the  algorithm  in  specific  situations  and  see  what  happens. 


The  steering  algorithm  GRAMPS  was  tested  for  several  dipole  distributions  with 
different  dipole  strengths  Kj  located  at  points  (xj,  yj,  Zi  =  0).  A  10  x  10  x  10  grid  was 

used.  The  missile  was  assumed  to  be  flying  at  a  predetermined  constant  velocity  vm  and  at 

=> 

time  t  =  0  to  be  located  at  point  (xQ  =  10,  y0  =  10,  Zq  =  10).  The  gradient  matrix  Q  was 
calculated  at  this  point  as  described  in  Section  II  with  GRAD  and  the  inverse  operation  was 
performed  with  code  GRAMPS.  One  of  the  two  resultant  vectors  pointing  in  the  forward 
hemisphere  was  selected  and  the  missile  moved  in  this  direction  for  time  period  At  to  point 
(x i ,  yj,  zi).  At  this  point  the  process  was  repeated  and  another  direction  selected 
corresponding  to  the  previous  vector  found  at  point  (x0,  y0,  z0).  This  process  was 
repeated  until  the  missile  either  hit  or  missed  the  target.  Table  2  gives  the  dipole  strengths 
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Table  2.  Dipole  Strength  and  Location  Used  for  the  Calculations  Shown 

In  Figs.  A-5  through  A-13 


Dipole  Location 

Results  Shown 

Dipole  Strengths 

(x,  y,  z) 

in  Figure  No. 

1.0 

(0,  0,  0) 

A-5 

1.0 

(0,  0,  0) 

A-6 

0.1 

(10,0,0) 

1.0 

(0,  0,  0) 

A-7 

0.5 

(10,  0,  0) 

1.0 

(0,  0,  0) 

A-8 

0.1 

(0,  10,  0) 

1.0 

(0,  0,  0) 

A-9 

0.2 

(0,  10,  0) 

1.0 

(0,  0,  0) 

A-10 

0.1 

(0,10, 0) 

0.1 

(5,  5,  0) 

1.0 

(0,  0,  0) 

A-1 1 

0.2 

(0,  10,  0) 

0.1 

(5,  5,  0) 

1.0 

(0,  0,  0) 

A-12 

0.1 

(0, 10,  0) 

0.2 

(5,  5,  0) 

1.0 

(0,0,  0) 

A-13 

0.2 

(0,  10,  0) 

0.2 

(5,  5,  0) 

and  locations  used  in  a  number  of  steering  calculations  and  Figs.  A-5  through  A- 15  show 
the  results  of  the  calculations.  In  each  figure,  drawing  (a)  shows  the  initial  vectors 

— ^  —4 

indicated  by  a,  b,  a',  and  b',  and  drawings  (b)  and  (c)  show  paths  taken  as  guided  by  the 
algorithm  GRAMPS.  Note  that  although  a  symmetric  situation  was  selected  setting  initial 
conditions  for  the  run  shown  in  Fig.  A-14,  the  algorithm  still  guides  the  missile  to  a  target. 
By  symmetry  consideration  one  might  expect  the  missile  to  be  guided  between  the  targets. 
It  seems,  however,  that  a  slight  numerical  error  breaks  the  symmetry  and  prevents  this 
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from  happening.  In  a  real  situation  in  the  field  we  expect  that  this  would  also  occur  and 
that  slight  perturbations  in  the  fields  would  break  the  symmetry  and  force  the  missile  to  one 
of  the  two  targets. 


Our  limited  sample  shows  that: 

1 .  The  steering  algorithm  works  when  the  target  is  isolated  from  decoys.  When 
decoys  are  present,  the  steering  algorithm  can  be  misled  as  shown  in  Fig.  10(c) 
and  (d).  In  this  case  the  decoy  of  strength  K2  =  10  was  10  units  away  from 


75 

the  starting  point  or  —  times  closer  than  the  target.  Since  the  dipole  field  goes 


like - 3  this  means  that  the  interaction  between  target  and  sensor  and 

(distance) 


72 

decoy  and  sensor  is  in  the  ratio  of  — .  The  decoy  has  a  stronger  attraction  for 
the  target. 

2.  Either  of  the  vectors  a  or  b  pointing  in  the  forward  hemisphere  will  direct  or 
steer  the  missile  to  the  desired  target. 


3 .  It  seems  that  the  size  of  the  angle  between  the  vectors  from  the  inverse  solution 
may  be  a  good  indication  of  the  steering  effectiveness.  If  the  two  vectors  that 
point  in  the  forward  direction  subtend  on  angle  <  90  deg  the  steering  will  be 
successful.  If  this  angle  is  >  90  deg  a  false  target  will  be  reached. 


Clearly,  further  analysis  is  required.  The  present  results  only  point  the  ways  to 
more  simulations  and  analysis.  These  results,  however,  do  indicate  feasibility  of  steering 
by  using  the  gradiometer  measurements. 
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V.  COUNTERMEASURES 


The  magnetic  SQUID  gradiometer  tank  seeker  could  be  overcome  by 

1 .  Building  tanks  out  of  a  nonmagnetic  alloy, 

2.  Reducing  the  signature  of  magnetized  tanks  (deperming  and  degaussing)  below 
a  measurable  level, 

3 .  Introducing  false  signatures, 

4.  Providing  decoys  to  simulate  the  tank’s  magnetic  signature  and  divert  the 
seeker. 

As  far  as  we  know  there  are  no  tanks  constructed  solely  from  nonmagnetic  material, 
and  it  is  not  likely  that  this  can  be  done  in  the  foreseeable  future.  Ceramic  armor  is  being 
used,  but  other  components  of  a  tank  are  not  built  completely  out  of  ceramic  or  non¬ 
magnetic  material. 

Degaussing  and  deperming  are  well  developed  procedures  used  by  the  Navy  to 
protect  surface  ships  and  submarines  from  mines.  Deperming  removes  some  of  the 
permanent  magnetization  of  the  vessel  produced  during  material  processing  and 
manufacture.  The  effect  lasts  for  about  two  years,  at  which  point  the  process  can  be 
repeated.  Degaussing  has  to  be  performed  continuously  because  it  removes  the  induced 
field  due  to  the  earth's  field.  This  has  to  be  done  on  demand  and  has  to  be  customized  to 
the  earth's  field  conditions. 

Degaussing  with  coils  on  a  surface  vessel  can  be  accomplished  easily  down  to 
20  percent  of  the  induced  magnetization  (vertical  component).  With  a  lot  more  effort,  a 
reduction  down  to  5  percent  is  predicted.  For  submarines,  degaussing  is  not  practiced 
because  of  the  weight  of  coils  and  power  requirements.  Deperming  is  more  practical  and  is 
used  for  both  surface  vessels  and  submarines. 

For  tanks,  no  operational  system  for  reducing  the  magnetic  signature  by  degaussing 
or  deperming  exists.  This  technology  has  been  evaluated  and  a  deperming  cage  constructed 
for  test  purposes.  Reduction  down  to  20  percent  of  original  value  has  been  obtained. 
Clearly,  this  number  could  be  improved,  but  the  cost  and  operational  utility  has  to  be 
established.  The  Army  seems  to  be  more  interested  in  signature  modification  designed  to 
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trigger  mines  well  ahead  of  tanks  than  in  signature  reduction.  A  lot  of  work  has  gone  into 
this  effort.  How  well  these  signature  modification  techniques  and  signal  reduction 
techniques  will  actually  work  in  the  field  to  countermeasure  the  seeker  has  to  be  determined 
through  tests,  evaluation,  and  further  analysis. 

Other  possible  countermeasures  include  the  use  of  permanent  or  electromagnets 
used  as  "flares"  to  distract  the  seeker.  For  the  permanent  magnet,  it  is  simple  to  compute 
the  required  size  of  the  decoy. 

Assume  that  the  decoy  is  made  from  a  magnetic  material  with  a  permanent 
magnetization  M  (B  =  47tM).  The  dipole  moment  of  such  a  decoy  is  independent  of  the 
details  of  the  shape  of  the  object  and  is  given  by  the  product  of  M  and  the  volume.  This 
can  be  seen  by  the  following  example.  Consider  a  cylinder  with  cross-sectional  area  A, 
and  length  L,  with  the  magnetization  M  oriented  along  the  cylinder  axis.  This  is  equivalent 
for  purposes  of  calculating  the  external  field,  to  a  surface  current  (in  cgs  units)  K  =  cM. 
The  dipole  moment  is  given  by  m  =  IA/c  =  KLA/c  =  MAL  =  M  Volume.  In  general,  the 
dipole  moment  of  the  material  is  given  by  m  =  j M  dV  =  M  Volume. 

The  dipole  moment  of  the  target  was  estimated  (conservatively)  as  m^get  = 
B0R3,  with  R  taken  as  1  m.  For  a  permanent  magnet  decoy,  m^ecoy  >  mtargei)  or  Mdccoy 
V decoy  >  B0R3.  For  Swedish  iron,  an  inexpensive  magnetic  material,  M  =  103 
corresponding  to  B  ~  12,000  gauss.  Since  BQR3is  approximately  unity,  V decoy  ~  10-3 
m3  =  103  cm3  (1  liter).  Each  of  these  decoys  would  therefore  weigh  about  8  kg. 
These  would  have  to  be  discharged  at  roughly  100-m  intervals  to  be  effective  as  decoys.  In 
10  km,  about  a  ton  of  decoys  would  be  needed.  More  exotic  material  may  raise  M  by  an 
order  of  magnitude  (decreasing  the  volume  by  the  same  amount),  but  would  increase  the 
cost. 

This  is  close  to  the  absolute  limit  to  be  expected  if  all  the  (relevant)  spins  in  the 
material  were  aligned.  The  magnetic  moment  per  spin  is  eh/47tmc  =  0.93  x  10~2°  gauss- 
cm3.  As  an  order  of  magnitude  estimate,  one  may  take  two  spins  per  atom  and  1023 
atoms/cm3,  this  yields  Miimit  =  1.8  x  103.  Therefore,  no  material  is  likely  to  extend  this 
limit  significantly.  The  use  of  an  electromagnet  is  unlikely  to  improve  the  values;  an  air 
core  magnet  will  have  a  lower  field  than  one  with  a  permeable  core,  and  these  will  be 
subject  to  the  same  limits  as  the  permanent  magnets. 


VI.  SUMMARY  AND  CONCLUSIONS 
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Ws  have  looked  at  the  feasibility  of  using  a  squid  gradiometer  as  a  sersor  for  an 
antitank  missile  terminal  guidance  system.  It  is  well  known  that  inversion  of  the 
gradiometer  signal  (from  a  collection  of  magnetic  dipoles)  gives  an  ambiguous  result  for  the 
direction  of  the  target;  four  vectors  in  pairs  pointing  in  opposite  directions.  For  this  reason 
the  determination  of  the  magnetic  gradient  in  general  would  not  be  a  very  useful 
measurement  for  a  tank  detection  and  guidance  system,  without  additional  information. 

In  our  application,  some  additional  information  is  available  and  two  directions  can 
be  selected  from  the  four  derived  vectors.  It  turns  out  that  both  vectors  lead  to  the  target  by 
different  routes  when  our  algorithm  is  applied.  We  have  performed  simulations  of  missile 
guidance  assuming  various  target  and  decoy  magnetic  strengths  and  spatial  distributions. 
How  well  the  algorithm  works  depends  on  the  relative  signal  strengths  and  distances  to 
missile  of  the  target  and  decoy.  Several  representative  situations  have  been  investigated. 

We  also  estimated  the  range  of  application  of  the  SQUID  gradiometer  system.  For 
an  isolated  tank  with  only  an  induced  moment  due  to  the  earth’s  field,  the  approximate 
distance  from  the  target  at  which  a  SQUID  gradiometer  can  begin  to  guide  the  missile  is 
150  m.  Thus,  for  the  guidance  system  to  be  effective,  the  missile  nas  to  be  shot  into  a 
"bucket"  300  m  in  diameter  with  the  target  in  the  center. 

Our  results  show  that  the  principle  of  inverting  a  gradiometer  signal  to  obtain 
guidance  for  a  missile  designed  to  attach  a  magnetic  target  is  sound.  Whether  the  system 
can  be  made  practical  depends  on  the  feasibility  of  building  a  SQUID  from  high 
temperature  superconductive  material  with  low  noise  characteristics  at  temperatures  high 
enough  so  that  lightweight  refrigeration  is  possible.  Since  work  on  SQUIDS  is  going  on 
now,  the  answers  to  these  questions  are  forthcoming. 

After  questions  regarding  SQUID  characteristics  and  refrigeration  requirements  are 
resolved,  questions  regarding  missile  dynamics  have  to  be  addressed.  Clearly,  the  missile 
speed  and  maneuverability  have  to  be  designed  around  the  sensor.  We  do  not  expect  that 


27 


the  sensor  could  be  merely  adapted  to  an  existing  missile.  Finally,  the  question 
countermeasures  has  been  addressed,  but  further  in-depth  work  in  this  area  has  to 
initiated. 
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1.  Distributed  Dipole  Field  Study.  Contour  plots  of  magnetic  field  components  [(as),  (b),  (c)]  and 
magnetic  field  gradient  components  [(d),  (e),  (f),  (g),  (h),  (i)],  plotted  in  the  x,  I,  z)  plane. 

Dipole  located  at  (0,0,0)  of  strength  1.0. 


re  A-2.  Distributed  Dipole  Field  Study.  Contour  plots  of  magnetic  field  components  [(as),  (b),  (c)]  and 
magnetic  field  gradient  components  [(d),  (e),  (f),  (g),  (h),  (I)],  plotted  In  the  x,  I,  2)  plcne. 

Dipoles  located  at  (0,0,0)  and  (0,10,0)  of  strengths  1.0  and  0.1,  respectively. 
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Igure  A-3.  Distributed  Dipole  Field  Study.  Contour  plots  of  magnetic  field  components  [(as),  (b),  (c)l  and 
magnetic  field  gradient  components  [(d),  (e),  (f),  (g),  (h),  (I)],  plotted  In  the  x,  I,  z)  plane. 

Dipoles  located  at  (0,0,0)  and  (10,0,0)  of  strengths  1.0  and  0.5,  respectively. 


4.  Distributed  Dipole  Field  Study.  Contour  plots  of  magnetic  field  components  [(as),  (b),  (c)]  and 
magnetic  field  gradient  components  [(d),  (e),  (f),  (g),  (h),  (l)],  plotted  In  the  x,  I,  z)  plane. 

Ipoles  located  at  (0,0,0),  (0,10,0)  and  (5,5,0)  of  strengths  1.0,  0.5,  and  0.1,  respectively. 


Figure  A-5.  Simulation  of  steering  by  gradlometer  Initial  conditions  with  derived 

vectors  i,  b,  a*  and  b’  are  shown  In  (a),  and  missile  trajectories  following  vectors 

a  and  b  are  shown  In  (b)  and  (c).  A  magnetic  dipole  of  strength  K  =  1.0 
Is  assumed  to  be  located  at  (0,0,0)  In  this  calculation. 


Figure  A-7.  Simulation  of  steering  by  gradlometer  Initial  conditions  with  derived 

•  vectors  a,  b,  a'  and  b’  are  shown  in  (a),  and  missile  trajectories  following  vectors 

a  and  b  are  shown  in  (b)  and  (c).  Magnetic  dipoles  of  strengths  Ki  =1.0  and 
K2  =  0.5  are  assumed  to  be  located  at  point  (0,0,0)  and  (10,0,0),  respectively. 
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Figure  A-8.  Simulation  of  steering  by  gradiometer  initial  conditions  with  derived 

vectors  a,  b,  a'  and  b'  are  shown  In  (a),  and  missile  trajectories  following  vectors 

a  and  b  are  shown  in  (b)  and  (c).  Magnetic  dipoles  of  strengths  K-t  =1.0  and 
K2  =  0.1  are  assumed  to  be  located  at  point  (0,0,0)  and  (0,10,0),  respectively. 


A-14 


Figure  A-9.  Simulation  of  steering  by  gradiometer  initial  conditions  with  derived 

vectors  a,  b,  a'  and  b’  are  shown  In  (a),  and  missile  trajectories  following  vectors 

a  and  b  are  shown  In  (b)  and  (c).  Magnetic  dipoles  of  strengths  Ki  =1.0  and 
Ko  =  -0.2  are  assumed  to  be  located  at  point  (0,0,0)  and  (0,10,0),  respectively. 


Figure  A-10.  Simulation  of  steering  by  gradiometer  Initial  conditions  with 
derived  vectors  a,  b,  a'  and  b'  are  shown  in  (a),  and  missile  trajectories 

following  vectors  a  and  b  are  shown  In  (b)  and  (c).  Magnetic  dipoles  of 
strengths  Ki  =  1.0,  K2  =  -0.1  and  K3  =  0.1  are  assumed  to  be  located 
at  point  (0,0,0),  (0,10,0)  and  (5,5,0),  respectively. 
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Figure  A-il.  Simulation  of  steering  by  gradlometer  initial  conditions  with 
derived  vectors  a,  b,  a'  and  b'  are  shown  in  (a),  and  missile  trajectories 

following  vectors  a  and  b  are  shown  In  (b)  and  (c).  Magnetic  dipoles  of 
strengths  K-i  =  1.0,  K2  s  -0.2  and  K3  =  0.1  are  assumed  to  be  located 
at  point  (0,0,0),  (0,10,0)  and  (5,5,0),  respectively. 
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Figure  A-12.  Simulation  of  steering  by  gradiometer  initial  conditions  with 

derived  vectors  a,  b,  a'  and  b'  are  shown  In  (a),  and  missile  trajectories  following 

.  — • 

vectors  a  and  b  are  shown  in  (b)  and  (c).  Magnetic  dipoles  of  strengths 
Ki  =  i.o,  K2  =  0.1  and  K3  =  0.2  are  assumed  to  be  located  at  point 
(0,0,0),  (0,10,0)  and  (5,5,0),  respectively. 


Figure  A-13.  Simulation  of  steering  by  gradiometer  Initial  conditions  with 
derived  vectors  a,  b,  a'  and  b’  are  shown  In  (a),  and  missile  trajectories  following 

vectors  a  and  b  are  shown  in  (b)  and  (c).  Magnetic  dipoles  of  strengths 
Ki  =  1.0,  K2  s  0.2  and  K3  =  0.2  are  assumed  to  be  located  at  point 
(0,0,0),  (0,10,0)  and  (5,5,0),  respectively. 


Figure  A-14.  Simulation  of  steering  by  gradiometer  initial  conditions  with 
derived  vectors  a,  b,  a'  and  b'  are  shown  In  (a),  and  missile  trajectories  following 

vectors  a  and  b  are  shown  in  (b)  and  (c).  Magnetic  dipoles  of  strengths 
Ki  =1.0  and  «2  =  1.0  are  assumed  to  be  located  at  point 
(0,-5,0)  and  (0,5,0),  respectively. 
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Figure  A-15.  Simulation  of  steering  by  gradiometer  Initial  conditions  with 

derived  vectors  a,  b,  a'  and  bf  are  shown  In  (a),  and  missile  trajectories  following 

_  — * 

vectors  a  and  b  are  shown  In  (b)  and  (c).  Magnetic  dipoles  of  strengths 
Ki  =  1.0  and  K2  =  0.2  are  assumed  to  be  located  at  point 
(0,-7, 0)  and  (0,3,0),  respectively. 


